Exact factorization of the time-dependent electron-nuclear wavefunction 



O 

(N 

m 



5-H 

CD 



I 

O 

o 



> 
oo 
cn 

O 
O 



X 



All Abedi/'^ Neepa T. Maitra,^ and E.K.U. Gross^'^ 

^ Max-Planck Institut fiir Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany 
^European Theoretical Spectroscopy Facility (ETSF) 
^Department of Physics and Astronomy, Hunter College and the City University of New York, 
695 Park Avenue, Neiv York, New York 10065, USA 
(Dated: September 24, 2010) 

We present an exact decomposition of the complete wavefunction for a system of nuclei and elec- 
trons evolving in a time-dependent external potential. We derive formally exact equations for the 
nuclear and electronic wavefunctions that lead to rigorous definitions of a time-dependent potential 
energy surface (TDPES) and a time-dependent geometric phase. For the molecular ion exposed 
to a laser field, the TDPES proves to be a useful interpretive tool to identify different mechanisms of 
dissociation. 
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Treating electron-ion correlations in molecules and 
solids in the presence of time-dependent external fields 
is a major challenge, especially beyond the perturba- 
tive regime. To make numerical calculations feasible, 
the description usually involves approximations such 
as classical dynamics for nuclei with electron-nuclear 
coupling provided by Ehrenfest dynamics or surface- 
hopping LU, or even just static nuclei [2]. Quantum 
features of the nuclear d5mamics (e.g., zero-point ener- 
gies, tunneling, and interference) are included approx- 
imately in some methods |3, 4], while numerically ex- 
act solutions of the time-dependent Schrodinger equa- 
tion (TDSE) for the coupled system of electrons and nu- 
clei have been given for very small systems like H2 [5] . 
Clearly, the full electron-nuclear wavefunction contains 
the complete information on the system, but it lacks the 
intuitive picture that potential energy surfaces (PES) can 
provide. To this end, approximate TDPES were intro- 
duced by Kono |@] as instantaneous eigenvalues of the 
electronic Hamiltonian, and proved extremely useful in 
interpreting system-field phenomena. The concept of a 
TDPES arises in a different way in Cederbaum's recent 
work, where the Born-Oppenheimer (BO) approxima- 
tion is generalized to the time-dependent case 

In the present Letter we provide a rigorous separa- 
tion of electronic and nuclear motion by introducing an 
exact factorization of the full electron-nuclear wavefimc- 
tion. The factorization is a natural extension of the work 
of Himter |@], in which an exact decomposition was de- 
veloped for the static problem. It leads to an exact def- 
inition of the TDPES as well as a^erry vector poten- 
tial. Berry-Pancharatnam phases |9] are usually inter- 
preted as arising from an approximate decoupling of a 
system from "the rest of the world", thereby making the 
system Hamiltonian dependent on some "environmen- 
tal" parameters. For example, in the static BO approx- 
imation, the electronic Hamiltonian depends paramet- 
rically on the nuclear positions; i.e., the stationary elec- 
tronic Schrodinger equation is solved for each fixed nu- 
clear configuration R, yielding R-dependent eigenval- 
ues (the BO PES) and eigenfunctions (the BO wavefunc- 



tions). If the total molecular wavefunction is approxi- 
mated by a single product of a BO wavefunction and 
a nuclear wavefunction, the equation of motion of the 
latter contains a Berry-type vector potential. One may 
ask: is the appearance of Berry phases a consequence 
of the BO approximation or does it survive in the ex- 
act treatment? In this Letter we demonstrate that even 
in the exact treatment of the electron-nuclear system a 
Berry connection appears and we prove a new relation 
between this connection and the nuclear velocity field. 
For a numerically exactly solvable system we calculate 
the exact TDPES, demonstrate their interpretive power, 
and compare with approximate treatments. Through- 
out this paper we use atomic units and the electronic 
and nuclear coordinates are collectively denoted by r, 
R. The Hamiltonian for a system of interacting electrons 
and nuclei, evolving under a time-dependent external 
potential, may be written as 



H = Hbo 



(1) 



where Hbo is the traditional BO electronic Hamiltonian, 

HBO=fe (r) + Wee (l) + Ven (l, E) + Wnn (E) • (2) 

Here Tn{Te) is the nuclear(electronic) kinetic energy 
operator, Wnn (Weei) is the nuclear-nuclear (electron- 
electron) interaction, and t^"j(R, t) and V^^^_{r,t) are 
time-dependent external potentials acting on the nu- 
clei and electrons, respectively. The complete electron- 
nuclear wavefunction satisfies the TDSE: 



i?*(r,R,t) = iat*(r,R,i). 



(3) 



The central statement of this Letter is the following: 

Theorem I. (a) The exact solution of Eq. ^ can be writ- 
ten as a single product 

*(i,E,t) = %(i,t)x(E,i) (4) 
where <I>R(r, t) satisfies the normalization condition. 
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(5) 



for any fixed nuclear configuration, R, at any time t. 
(b) The wavefunctions $R,(r, t) and x(R, t) satisfy: 

[Hei{L,R,t) - e®,t))%(l,t) = idt<P^ig,t), (6) 

(E^HV. + A.(R,t)f 

- (7) 

where the electronic Hamiltonian is 



Hel (l, E, =ffBO + Kxt (l, + 5Z ]^ 



rHV,- A,(E,t))' , /-iV,x (8) 



2 V X 

+ A^(E, (-*V^ - A^(E, t)) 

Here the scalar and vector potential terms are 

e(R,t) = Unit) Hei{{r,R,t) - idt $r(<)\ (9) 



A.(R,t) = ($R(t) 



(10) 



where denotes an inner product over all elec- 

tronic variables only. 

Proof. Part (a): We must show that the exact solu- 
tion vl/(r, R, t) of the full TDSE (|3]l can be factorized as 
in Eqs. (|4Jl-(|5ll. To show this, choose at each point in 

time x(E, t) = e''^®'*) y^/ i)P and $r(i, t) = 

E, t) /x(E, t)i where ^(E, t) is real. The normaliza- 
tion condition (|5]l then follows immediately. 

Part (b): To derive Eqs. ll6|l- ((T0l l, we apply Frenkel's 
stationary action principle, 5 J^^ dt{^\H — idt\^) = , 
to the wavef unction @. We require the action to be sta- 
tionary with respect to variations in $R(r, t) and x(Ei t)f 
subject to the condition (|5]l. This then leads, after some 
algebra, to Eqs. l(6ll- ((T0l l. Hence, the product wavefunc- 
tion (ID is a stationary point of the action functional, but 
we still have to prove that this stationary point corre- 
sponds to an exact solution of the TDSE. By evaluating 
i9t($R(r, i)x(E5 1)) arid inserting Eqs. ©-([Zll, we verify 
that theluU TDSE (3) is satisfied. 

Theorem II. (a) Eqs. l|6]l-l|8ll are form-invariant under the 
following gaugelike transformation 

^ = exp(i0(E, <))%(!, 

x(E, t) ^ x(E, t) = cxp(-z0(E t))x(E, t) (11) 



A^(E,t) A^(E,i) = A^(E,i) + V^6'(E,i) 
e(E, t) e(E, t) = e(E, t) + dtOiR, t) (12) 



(b) The wavef imctions <I>R(r, t) and x(Eii) yielding 
a given solution, (r, Ei^)/ of Eq. lO are unique up to 
within the (Ei t)-dependent phase transformation ifTTI l. 

Proof. The form invariance of Eqs. l|6]l-l|8| is easily ver- 
ified by inserting Eqs. l(TT1l - l(T2l l into Eqs. (|6]l-||8| which 
proves part (a). To prove part (b), assume the exact 
wavefunction can be represented by two different prod- 
ucts: *(i,Ei) = %(l,i)x(E,i) = %(l,Ox(E,i)- 
Defining ,9(E,i) := x(E, 0/x(E, 0. then = 
l.9(E7 OPI*J'r(Ii Integrating this over r and using 
Eq. ©, we get |.g(E,OP = 1 implying g(E,0 = e*^®'*^ 
and hence the desired result $R(r, t) = e*^'=*'$R(r, t). 

The wavefunctions <l'R(r, t) and xCEi have a clear- 
cut physical meaning: |x(E)i)P ~ / l*(liEii)Pc^r is 
the probability density of finding the nuclear configura- 
tion E at time t, and |$r(£, = |«'(i,E,<)lVlx(E,^)P 
is the conditional probability of finding the electrons at 
r, given that the nuclear configuration is E- At locations 
where IxlEi approaches zero the TDPES may show 
peaks, in close analogy to the "quantum potential" in 
the Bohmian formulation of quantum mechanics. 
Eqs. ll6ll- (|T0l l determine the exact time-dependent molec- 
ular wavefunction, given an initial state. As writ- 
ten, the nuclear equation is particularly appealing as 
a Schrodinger equation with both scalar and vector- 
potential coupling terms contributing effective forces 
on the nuclei including any geometric phase effects. 
We call e(E,i) and A(E,i) the exact TDPES and time- 
dependent Berry connection, respectively. These two 
quantities mediate the coupling between the nuclear 
and the electronic degrees of freedom in a formally ex- 
act way. Eqs. |[6ll-([T0ll demonstrate that a Berry connec- 
tion indeed appears in the exact treatment. But does it 
produce a real effect or can it be gauged away by a suit- 
able choice of f?(E, t) in Eqs. ((TT|l - ((T2)l ? To shed some 
light on this question, we now prove an alternate ex- 
pression for the vector potential. Inserting $r = ^'/x 
into Eq. (lOt , and evaluating the nuclear gradient on this 
quotient, reveals that it is the difference of paramagnetic 
nuclear velocity fields derived from the full and nuclear 
wavefunctions : 

A.(E,t) = ^^-r^ ^- ,:r - (13) 



lx(E,i)l= 



lx(E,t)l= 



This equation is interesting in several respects. First, 

writing x(E,^) = e''^^^'*^| x(E , ^)|/ the last term on 
the right-hand-side of Eq. (|131 can be represented as 
Viy5(E, i)/ so it can be gauged away. Consequently, 
any true Berry connection (that cannot be gauged away) 
must come from the first term. If the exact ^'(t) is real- 
valued (e.g. for a non-current-carrying groimd state) 
then the first term on the right-hand-side of Eq. l(l3l l 
vanishes and hence the exact Berry connection vanishes. 
Second, since Im {^{t)\V ^'^{t))^ is the true nuclear 

(many-body) current density, Eq. (|T3] | implies that the 
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gauge-invariant current density, Im{x*'^ vX) + Ixl^A^, 
that follows from Eq. ((Tjl does indeed reproduce the ex- 
act nuclear current density (To|] . Hence, the solution 
x(R, t) of Eq. ^ is, in every respect, the proper nuclear 
many-body wavefunction: Its absolute-value squared 
gives the exact nuclear (A^-body) density while its phase 
yields the correct nuclear (A^-body) current density. 

In the following, we first discuss some limiting 
cases of the exact Eqs. |[6ll-<[T0ll. Fixing the gauge via 

^R(i) [^{^^.(t) ) = 0, the electronic equation reads 



(14) 



<I>R while the nuclear 
Eq~0 with e(E,0 = 
Note that the electronic 



with (/)R = 
equation retains its form 

Eq. l(T4l l and the nuclear Eq. ||7ll have to be prop- 
agated simultaneously because the Hamiltonian (|8) 
depends on x®!^)- Taking the large-nuclear-mass 
limit of Eq. (|8), the electronic Hamiltonian reduces to 
Hei > Hbo + y^xt' the dependence of i^ei on 

X drops out and the electronic Eq. |[T4t depends on the 
nuclear configuration R only parametrically. So, in this 
limit, Eq. ifMl l is propagated in time for each fixed nu- 
clear configuration R, which is precisely Cederbaum's 
time-dependent generalization of the BO approxima- 
tion 0]. Hence, the full Eqs. ©-lO can be viewed as an 
"exactification" of the intuitively appealing procedure 
of Ref. 10] ■ If/ furthermore, we treat the nuclei classi- 
cally, i.e., use the Hamiltonian (0 to generate classical 
equations of motion for the nuclei, we obtain 



MR 



E„ + R X B„ 



(15) 



where the electric and magnetic "Berry fields" are given 
by = V^e(E,t) - ^ and = x A^. The 
additional magnetic field was also found, in the appro- 
priate limit, in an exact path-integral approach to the 
coupled dynamics fllll , and also in other work fT2ll . Be- 
ing strictly equivalent to the TDSE, the electronic and 
nuclear Eqs. ©-lIZl provide a rigorous starting point 
suitable for making systematic semiclassical approxima- 
tions jl^[3 beyond the purely classical limit of Eq. lITSl l. 

We now return to the exact formulation to investigate 
the TDPES for a numerically exactly solvable model: the 
H2 molecular ion subject to a linearly polarized laser 
field. By restricting the motion of the nuclei and the elec- 
tron to the direction of the polarization axis of the laser 
field , the problem can be modelled with a ID Hamilto- 
nian featuring "soft-Coulomb" interactions | [l5l| : 



H{t) 



1 32 



1 32 



1 



M ai?2 2^i, dz^ V0.03 + i?2 
1 1 



+ Vi{z,t) 
(16) 
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FIG. 1. Snapshots of the TDPES (blue lines) and nuclear den- 
sity (black) at times indicated, for the molecule subject to 
the laser-field (see text), h — 10^''W/cm^ (dashed line) and 
I2 = 2.5 X (solid line). The circles indicate the po- 

sition and energy of the classical particle in the exact-Ehrenfest 
calculation (Ji: open, J2: solid). For reference, the ground-state 
BO surface is shown as the thin red line. 



where R and z are the rnternuclear distance and the elec- 
tronic coordinate as measured from the nuclear center- 
of-mass, respectively, and the electronic reduced mass is 
given by /ig = (2i\/)/(2M+l),M being the proton mass. 
The laser field is represented by Vi{z, t) = qezE{t) where 
E{t) denotes the electric field amplitude and the re- 
duced charge Qe = (2M+2) /{2M +1). We consider a A = 
228 nm laser field, represented by E{t) = Eof{t) sm{Lut), 
for two peak intensities, Ii ~ |£'oP = 10^^VK/cm2 and 
h = l^'oP = 2.5 X 10^^VF/cm2. The envelope function 
f{t) is chosen such that the field is linearly ramped from 
zero to its maximum strength at ^ = 7.6 fs and thereafter 
held constant. 

Starting from the exact ground-state as initial con- 
dition, we propagate the TDSE numerically to obtain 
the full molecular wavefunction \I'(z,i?,t). As there 
is only one nuclear degree of freedom (after sepa- 
rating off the center-of-mass motion), we can fix the 
gauge in Eqs. ((TTl l- l(T2t such that the vector poten- 
tial ( [T3t is always zero. From the computed exact 
time-dependent molecular wavefunction we compute 
the TDPES's; these, along with the corresponding nu- 
clear density, \x{R, are plotted in Fig.[T]at six snap- 
shots of time. The initial TDPES lies practically on top 
of the ground-state BO surface, which is plotted in all 
the snapshots for comparison. Fig. |2] shows the exact 

internuclear distance (^'(t)| R I'^it)), along with the re- 
sults from three approximate methods: (i) the usual 
Ehrenfest approximation (i.e. Eq. [TSl l, (ii) the "exact- 
Ehrenfest" approximation, which substitutes the exact 
TDPES for the Ehrenfest potential in the usual Ehren- 
fest approach and, (iii) an uncorrelated approach, the 
time-dependent Hartree (self-consistent field) approx- 
imation, 4'ff(r, R, = 0(r, i)x(Ri 0' where the elec- 
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tronic part does not depend on R at all. Fig.|2]shows that 
for the intensity Ii, all methods yield dissociation, while 
for the weaker I2, only the exact does. We now discuss 
how the TDPES contains the signature of this behavior. 
Note that the laser-field does not couple directly to the 
nuclear relative coordinate R, but only indirectly via the 
TDPES. 

= \Eo\'^ = IQ^^W/cm^: The dissociation of the 
molecule is dramatically reflected in the exact TDPES, 
whose well flattens out, causing the nuclear density to 
spill to larger separations. Importantly, the tail of the 
TDPES alternately falls sharply and returns in corre- 
spondence with the field, letting the density out; the TD- 
PES is the only potential acting on the nuclear system 
and transfers energy from the accelerated electron to the 
nuclei. The expectation value of the internuclear dis- 
tance in Fig. |2l demonstrates that among all the approx- 
imate calculations employed here, the exact-Ehrenfest is 
most accurate. Surprisingly, it even does better than TD- 
Hartree which treats the protons quantum mechanically, 
thus showing the importance of electron-nuclear corre- 
lation. 

h = \Eo\^ = 2.5 X lO^^W/cm^ : From Fig.El the exact 
calculation leads to dissociation, while none of the ap- 
proximations do, in contrast to the previous case. The 
TDPES of Fig. [H suggests that tunneling is the lead- 
ing mechanism for the dissociation: a well remains at 
all times that traps a classical particle, which would os- 
cillate inside it, as indeed reflected in Fig. |21 (See also 
the solid circles in Fig.[lJ. Although the tail has similar 
oscillations as for Ii, this does not lead to dissociation 
of classical nuclei due to the barrier; the TDPES in this 
case transfers the field energy to the nuclei via tunnel- 
ing. Although the exact-Ehrenfest has a larger ampli- 
tude of oscillation than the others, it ultimately cannot 
tunnel through the barrier. 

This example demonstrates how studying the TD- 
PES reveals the mechanism of dissociation. Because the 
TDPES includes the electron-nuclear correlation exactly, 
we believe the exact-Ehrenfest dynamics is the best one 
could do within a classical treatment of the nuclei. There 
is a need to go beyond classical djmamics when the dis- 



sociation proceeds mainly via turmeling. 
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FIG. 2. The internuclear separation {R){t) for the same inten- 
sities as in Fig.l. Left panel : h. Right panel : h- 

In conclusion, we have presented a rigorous factor- 
ization of the complete molecular wavefunction into 
an electronic contribution, <&R(r, t), and a nuclear part, 
x(R, t). The exact nuclear A^-body density is |x(B:7 0P 
while |$R(r, t)p represents the conditional probability 
of finding the electrons at r, given the nuclear configu- 
ration R. Their exact equations of motion are deduced. 
Via these equations, the time-dependent potential en- 
ergy surface (|9]l and the time-dependent Berry connec- 
tion (|T0] | are defined as rigorous concepts. We demon- 
strated with numerical examples that the TDPES is a 
powerful tool to analyze and interpret different types 
of dissociation processes (direct vs tunneling). The ex- 
act splitting of electronic and nuclear degrees of free- 
dom presented by Eqs. (|6|l- l(T0l l lends itself as a rigorous 
starting point for making approximations, especially for 
the systematic development of semiclassical approxima- 
tions. As a first step we have shown how the Ehrenfest 
equations with Berry potential emerge from treating the 
nuclei classically in the large-nuclear-mass limit. 
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